Picosecond quantum-classical dynamics reveals that the coexistence of light-induced microbial and animal chromophore rotary motion modulates the isomerization quantum yield of heliorhodopsin.

Rhodopsins are light-responsive proteins forming two vast and evolutionary distinct superfamilies whose functions are invariably triggered by the photoisomerization of a single retinal chromophore. In 2018 a third widespread superfamily of rhodopsins called heliorhodopsins was discovered using functional metagenomics. Heliorhodopsins, with their markedly different structural features with respect to the animal and microbial superfamilies, offer an opportunity to study how evolution has manipulated the chromophore photoisomerization to achieve adaptation. One question is related to the mechanism of such a reaction and how it differs from that of animal and microbial rhodopsins. To address this question, we use hundreds of quantum-classical trajectories to simulate the spectroscopically documented picosecond light-induced dynamics of a heliorhodopsin from the archaea thermoplasmatales archaeon (TaHeR). We show that, consistently with the observations, the trajectories reveal two excited state decay channels. However, inconsistently with previous hypotheses, only one channel is associated with the -C13C14- rotation of microbial rhodopsins while the second channel is characterized by the -C11C12- rotation typical of animal rhodopsins. The fact that such -C11C12- rotation is aborted upon decay and ground state relaxation, explains why illumination of TaHeR only produces the 13-cis isomer with a low quantum efficiency. We argue that the documented lack of regioselectivity in double-bond excited state twisting motion is the result of an "adaptation" that could be completely lost via specific residue substitutions modulating the steric hindrance experienced along the isomerization motion.

Ultrafast excited-state dynamics of Luteins in the major light-harvesting complex LHCII.

Carotenoid pigments are known to present a functional versatility when bound to light-harvesting complexes. This versatility originates from a strong correlation between a complex electronic structure and a flexible geometry that is easily tunable by the surrounding protein environment. Here, we investigated how the different L1 and L2 sites of the major trimeric light-harvesting complex (LHCII) of green plants tune the electronic structure of the two embedded luteins, and how this reflects on their ultrafast dynamics upon excitation. By combining molecular dynamics and quantum mechanics/molecular mechanics calculations, we found that the two luteins feature a different conformation around the second dihedral angle in the lumenal side. The s-cis preference of the lutein in site L2 allows for a more planar geometry of the π -conjugated backbone, which results in an increased degree of delocalization and a reduced excitation energy, explaining the experimentally observed red shift. Despite these remarkable differences, according to surface hopping simulations the two luteins present analogous ultrafast dynamics upon excitation: the bright S 2 state quickly decays (in ∼ 50 fs) to the dark intermediate S x , eventually ending up in the S 1 state. Furthermore, by employing two different theoretical approaches (i.e., Förster theory and an excitonic version of surface hopping), we investigated the experimentally debated energy transfer between the two luteins. With both approaches, no evident energy transfer was observed in the ultrafast timescale.

Triterpene Derivatives as Potential Inhibitors of the RBD Spike Protein from SARS-CoV-2: An In Silico Approach.

The appearance of a new coronavirus, SARS-CoV-2, in 2019 kicked off an international public health emergency. Although rapid progress in vaccination has reduced the number of deaths, the development of alternative treatments to overcome the disease is still necessary. It is known that the infection begins with the interaction of the spike glycoprotein (at the virus surface) and the angiotensin-converting enzyme 2 cell receptor (ACE2). Therefore, a straightforward solution for promoting virus inhibition seems to be the search for molecules capable of abolishing such attachment. In this work, we tested 18 triterpene derivatives as potential inhibitors of SARS-CoV-2 against the receptor-binding domain (RBD) of the spike protein by means of molecular docking and molecular dynamics simulations, modeling the RBD S1 subunit from the X-ray structure of the RBD-ACE2 complex (PDB ID: 6M0J). Molecular docking revealed that at least three triterpene derivatives of each type (i.e., oleanolic, moronic and ursolic) present similar interaction energies as the reference molecule, i.e., glycyrrhizic acid. Molecular dynamics suggest that two compounds from oleanolic and ursolic acid, OA5 and UA2, can induce conformational changes capable of disrupting the RBD-ACE2 interaction. Finally, physicochemical and pharmacokinetic properties simulations revealed favorable biological activity as antivirals.

How the pH Controls Photoprotection in the Light-Harvesting Complex of Mosses.

In response to varying light conditions, light-harvesting complexes (LHCs) switch from a light-harvesting state to a quenched state to protect the photosynthetic organism from excessive light irradiation in a strategy known as nonphotochemical quenching (NPQ). NPQ is activated by an acidification of the thylakoid lumen, which is sensed directly or indirectly by the LHC, resulting in a conformational change of the complex that leads to the quenched state. The conformational changes responsible for NPQ activation and their connection to specific quenching mechanisms are still unknown. Here, we investigate the pH-triggered conformational changes in the light-harvesting complex stress-related (LHCSR) of mosses. By combining constant-pH molecular dynamics and enhanced sampling techniques, we find that the pH sensitivity of the complex is driven by the coupled protonation of three residues modulating the conformation of the short amphipathic helix placed at the lumen side of the embedding membrane. Combining these results with quantum mechanics/molecular mechanics calculations, we show that the quenching mechanism sensitive to the pH goes through a charge-transfer between a carotenoid and an excited chlorophyll, which is controlled by the protein conformation.

Automated QM/MM Screening of Rhodopsin Variants with Enhanced Fluorescence.

We present a computational protocol for the fast and automated screening of excited-state hybrid quantum mechanics/molecular mechanics (QM/MM) models of rhodopsins to be used as fluorescent probes based on the automatic rhodopsin modeling protocol (a-ARM). Such "a-ARM fluorescence screening protocol" is implemented through a general Python-based driver, PyARM, that is also proposed here. The implementation and performance of the protocol are benchmarked using different sets of rhodopsin variants whose absorption and, more relevantly, emission spectra have been experimentally measured. We show that, despite important limitations that make unsafe to use it as a black-box tool, the protocol reproduces the observed trends in fluorescence and it is capable of selecting novel potentially fluorescent rhodopsins. We also show that the protocol can be used in mechanistic investigations to discern fluorescence enhancement effects associated with a near degeneracy of the S1/S2 states or, alternatively, with a barrier generated via coupling of the S0/S1 wave functions.

The three-center two-positron bond.

Computational studies have shown that one or more positrons can stabilize two repelling atomic anions through the formation of two-center positronic bonds. In the present work, we study the energetic stability of a system containing two positrons and three hydride anions, namely 2e+[H3 3-]. To this aim, we performed a preliminary scan of the potential energy surface of the system with both electrons and positrons in a spin singlet state, with a multi-component MP2 method, that was further refined with variational and diffusion Monte Carlo calculations, and confirmed an equilibrium geometry with D 3h symmetry. The local stability of 2e+[H3 3-] is demonstrated by analyzing the vertical detachment and adiabatic energy dissociation channels. Bonding properties of the positronic compound, such as the equilibrium interatomic distances, force constants, dissociation energies, and bonding densities are compared with those of the purely electronic H3 + and Li3 + systems. Through this analysis, we find compelling similarities between the 2e+[H3 3-] compound and the trilithium cation. Our results strongly point out the formation of a non-electronic three-center two-positron bond, analogous to the well-known three-center two-electron counterparts, which is fundamentally distinct from the two-center two-positron bond [D. Bressanini, J. Chem. Phys., 2021, 155, 054306], thus extending the concept of positron bonded molecules.

Retinal chromophore charge delocalization and confinement explain the extreme photophysics of Neorhodopsin.

The understanding of how the rhodopsin sequence can be modified to exactly modulate the spectroscopic properties of its retinal chromophore, is a prerequisite for the rational design of more effective optogenetic tools. One key problem is that of establishing the rules to be satisfied for achieving highly fluorescent rhodopsins with a near infrared absorption. In the present paper we use multi-configurational quantum chemistry to construct a computer model of a recently discovered natural rhodopsin, Neorhodopsin, displaying exactly such properties. We show that the model, that successfully replicates the relevant experimental observables, unveils a geometrical and electronic structure of the chromophore featuring a highly diffuse charge distribution along its conjugated chain. The same model reveals that a charge confinement process occurring along the chromophore excited state isomerization coordinate, is the primary cause of the observed fluorescence enhancement.

On the fluorescence enhancement of arch neuronal optogenetic reporters.

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter.

Evolution of the Automatic Rhodopsin Modeling (ARM) Protocol.

In recent years, photoactive proteins such as rhodopsins have become a common target for cutting-edge research in the field of optogenetics. Alongside wet-lab research, computational methods are also developing rapidly to provide the necessary tools to analyze and rationalize experimental results and, most of all, drive the design of novel systems. The Automatic Rhodopsin Modeling (ARM) protocol is focused on providing exactly the necessary computational tools to study rhodopsins, those being either natural or resulting from mutations. The code has evolved along the years to finally provide results that are reproducible by any user, accurate and reliable so as to replicate experimental trends. Furthermore, the code is efficient in terms of necessary computing resources and time, and scalable in terms of both number of concurrent calculations as well as features. In this review, we will show how the code underlying ARM achieved each of these properties.

Pro219 is an electrostatic color determinant in the light-driven sodium pump KR2.

Color tuning in animal and microbial rhodopsins has attracted the interest of many researchers, as the color of their common retinal chromophores is modulated by the amino acid residues forming the chromophore cavity. Critical cavity amino acid residues are often called "color switches", as the rhodopsin color is effectively tuned through their substitution. Well-known color switches are the L/Q and A/TS switches located in the C and G helices of the microbial rhodopsin structure respectively. Recently, we reported on a third G/P switch located in the F helix of the light-driven sodium pumps of KR2 and JsNaR causing substantial spectral red-shifts in the latter with respect to the former. In order to investigate the molecular-level mechanism driving such switching function, here we present an exhaustive mutation, spectroscopic and computational investigation of the P219X mutant set of KR2. To do so, we study the changes in the absorption band of the 19 possible mutants and construct, semi-automatically, the corresponding hybrid quantum mechanics/molecular mechanics models. We found that the P219X feature a red-shifted light absorption with the only exception of P219R. The analysis of the corresponding models indicate that the G/P switch induces red-shifting variations via electrostatic interactions, while replacement-induced chromophore geometrical (steric) distortions play a minor role. However, the same analysis indicates that the P219R blue-shifted variant has a more complex origin involving both electrostatic and steric changes accompanied by protonation state and hydrogen bond networks modifications. These results make it difficult to extract simple rules or formulate theories for predicting how a switch operates without considering the atomistic details and environmental consequences of the side chain replacement.

Frontiers in Multiscale Modeling of Photoreceptor Proteins.

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Covalent bonds in positron dihalides.

We report a computational study on homo- and heteronuclear e+[X-Y-] compounds formed by two halide anions (X-, Y- = F-, Cl-, Br-) and one positron. Our results indicate the formation of energetically stable positronic molecules in all cases. Analysis of the electron and positron densities points out that the formation of positron covalent bonds underlies the stabilization of the otherwise repelling dihalides, revealing that positronic bonding can reach far beyond the previously addressed e+[H-H-] molecule [J. Charry, M. T. do N. Varella and A. Reyes, Angew. Chem. Int. Ed., 2018, 57, 8859-8864.]. To a significant extent, the properties of the positron dihalides are similar to those of the purely electronic analogs, e-[A+B+], molecular cations with isoelectronic atomic cores (A+, B+ = Na+, K+, Rb+) bound by one electron. The positron bonds in the e+[X-Y-] complexes are however stronger than those in the isoelectronic e-[A+B+] counterparts, as the former have shorter bond lengths and higher bond energies. While an energy decomposition analysis points out that both electronic and positronic bonds essentially arise from electrostatic interactions, the more stable positron bonds are partly due to the higher polarizabilities of the dihalide anions, and partly to more significant contributions from correlation and relaxation effects.

Modern quantum chemistry with [Open]Molcas.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Web-ARM: A Web-Based Interface for the Automatic Construction of QM/MM Models of Rhodopsins.

This article introduces Web-ARM, a specialized tool, online available, designed to build quantum mechanical/molecular mechanical models of rhodopsins, a widely spread family of light-responsive proteins. Web-ARM allows the rapidly building of models of rhodopsins with a documented quality and the prediction of trends in UV-vis absorption maximum wavelengths, based on their excitation energies computed at the CASPT2//CASSCF/Amber level of theory. Web-ARM builds upon the recently reported, python-based a-ARM protocol [J. Chem. Theory Comput., 2019, 15, 3134-3152] and, as such, necessitates only a crystallographic structure or a comparative model in PDB format and a very basic knowledge of the studied rhodopsin system. The user-friendly web interface uses such input to generate congruous, gas-phase models of rhodopsins and, if requested, their mutants. We present two possible applications of Web-ARM, which showcase how the interface can be employed to assist both research and educational activities in fields at the interface between chemistry and biology. The first application shows how, through Web-ARM, research projects (e.g., rhodopsin and rhodopsin mutant screening) can be carried out in significantly less time with respect to using the required computational photochemistry tools via a command line. The second application documents the use of Web-ARM in a real-life educational/training activity, through a hands-on experience illustrating the concepts of rhodopsin color tuning.

Multistate Multiconfiguration Quantum Chemical Computation of the Two-Photon Absorption Spectra of Bovine Rhodopsin.

Recently, progress in IR sources has led to the discovery that humans can detect infrared (IR) light. This is hypothesized to be due to the two-photon absorption (TPA) events promoting the retina dim-light rod photoreceptor rhodopsin to the same excited state populated via one-photon absorption (OPA). Here, we combine quantum mechanics/molecular mechanics and extended multiconfiguration quasi-degenerate perturbation theory calculations to simulate the TPA spectrum of bovine rhodopsin (Rh) as a model for the human photoreceptor. The results show that the TPA spectrum of Rh has an intense S0 → S1 band but shows also S0 → S2 and S0 → S3 transitions whose intensities, relative to the S0 → S1 band, are significantly increased when compared to the corresponding bands of the OPA spectrum. In conclusion, we show that IR light in the 950 nm region can be perceived by rod photoreceptors, thus supporting the two-photon origin of the IR perception. We also found that the same photoreceptor can perceive red (i.e., close to 680 nm) light provided that TPA induces population of S2.

a-ARM: Automatic Rhodopsin Modeling with Chromophore Cavity Generation, Ionization State Selection, and External Counterion Placement.

The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol-1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol-1 (0.04 eV).

Fast and accurate prediction of proton affinities: revisiting the extended Koopmans' theorem for protons.

In this work we propose schemes based on the extended Koopmans' theorem for quantum nuclei (eKT), in the framework of the any particle molecular orbital approach (APMO/KT), for the quantitative prediction of gas phase proton affinities (PAs). The performance of these schemes has been tested on a set of 300 organic molecules containing diverse functional groups. The APMO/KT scheme scaled by functional group (APMO/KT-SC-FG) displays an overall mean absolute error of 1.1 kcal mol-1 with respect to experimental data. Its performance in PA calculations is similar to that of post-Hartree-Fock composite methods or that of the APMO second order proton propagator (APMO/PP2) approach. The APMO/KT-SC-FG scheme is also employed to predict PAs of polyfunctional molecules such as the Nerve Agent VX and the 20 common α-amino acids, finding excellent agreement with available theoretical and/or experimental data. The accuracy of the predictions demonstrates that the APMO/KT-SC-FG scheme is a low-cost alternative to adiabatic methods for the calculation of accurate PAs. One of the most appealing features of the APMO/KT-SC-FG scheme, is that PAs can be derived from one single-point APMO Hartree-Fock calculation.

On the physical interpretation of the nuclear molecular orbital energy.

Recently, several groups have extended and implemented molecular orbital (MO) schemes to simultaneously obtain wave functions for electrons and selected nuclei. Many of these schemes employ an extended Hartree-Fock approach as a first step to find approximate electron-nuclear wave functions and energies. Numerous studies conducted with these extended MO methodologies have explored various effects of quantum nuclei on physical and chemical properties. However, to the best of our knowledge no physical interpretation has been assigned to the nuclear molecular orbital energy (NMOE) resulting after solving extended Hartree-Fock equations. This study confirms that the NMOE is directly related to the molecular electrostatic potential at the position of the nucleus.

Prediction of proton affinities of organic molecules using the any-particle molecular-orbital second-order proton propagator approach.

We assess the performance of the recently developed any-particle molecular-orbital second-order proton propagator (APMO/PP2) scheme [M. Díaz-Tinoco, J. Romero, J. V. Ortiz, A. Reyes and R. Flores-Moreno, J. Chem. Phys., 2013, 138, 194108] on the calculation of gas phase proton affinities (PAs) of a set of 150 organic molecules comprising several functional groups: amines, alcohols, aldehydes, amides, ketones, esters, ethers, carboxylic acids and carboxylate anions. APMO/PP2 PAs display an overall mean absolute error of 0.68 kcal mol-1 with respect to experimental data. These results suggest that the APMO/PP2 method is an alternative approach for the quantitative prediction of gas phase proton affinities. One novel feature of the method is that a PA can be obtained from a single calculation of the optimized protonated molecule.

Revisiting the Dielectric Constant Effect on the Nucleophile and Leaving Group of Prototypical Backside SN2 Reactions: A Reaction Force and Atomic Contribution Analysis.

The solvent effect on the nucleophile and leaving group atoms of the prototypical F- + CH3Cl → CH3F + Cl- backside bimolecular nucleophilic substitution reaction (SN2) is analyzed employing the reaction force and the atomic contributions methods on the intrinsic reaction coordinate (IRC). Solvent effects were accounted for using the polarizable continuum solvent model. Calculations were performed employing 11 dielectric constants, ε, ranging from 1.0 to 78.5, to cover a wide spectrum of solvents. The reaction force data reveal that the solvent mainly influences the region of the IRC preceding the energy barrier, where the structural rearrangement to reach the transition state occurs. A detailed analysis of the atomic role in the reaction as a function of ε reveals that the nucleophile and the carbon atom are the ones that contribute the most to the energy barrier. In addition, we investigated the effect of the choice of nucleophile and leaving group on the ΔE0 and ΔE‡ of Y- + CH3X → YCH3 + X- (X, Y = F, Cl, Br, I) in aqueous solution. Our analysis allowed us to find relationships between the atomic contributions to the activation energy and leaving group ability and nucleophilicity.